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Abstract 



Cosmological magnetic fields in open Friedmann universes can experience superadiabatic 
' amplification within the realm of conventional electromagnetism. This is possible mathemati- 

cally, despite the conformal invariance of Maxwell's equations, because Friedmann spacetimes 
with non-Euclidean spatial geometry are not globally conformal to Minkowski space. Physi- 
cally, this means that even universes that are marginally open today can sustain large-scale 
magnetic fields that are substantially stronger than previously anticipated. In the present 
article, we investigate this purely geometric amplification mechanism in greater detail, fo- 
cusing on the early evolution of the electromagnetic modes in inflationary Friedmann models 
with hyperbolic spatial geometry. This also allows us to refine the earlier numerical estimates 
1^ ■ and provide the current spectrum of the residual, superadiabatically amplified magnetic field. 

^ . PACS numbers: 98.80.-k, 98.62.En, 98.65.Dx 



1 Introduction 

Large-scale magnetic (B) fields have been observed throughout the universe: from galaxies and 
galaxy clusters to remote protogalactic structures pp. Nevertheless, the origin of these fields is 
still an open question. The galactic dynamo can explain the galactic and (possibly) the cluster 
B-fields, but the mechanism faces difficulties with those seen in high-redshift protogalaxies [2]. 
The problem deepens when one takes into account recent claims for the detection of magnetic 
fields in empty intergalactic space [3], where the dynamo amplification presumably cannot op- 
erate. Cosmology could provide the answer, but generating sustainable -B-fields in the early 
universe has been proved very challenging [JJ. The main obstacle is that conventional magnetic 
fields in spatially flat Friedmann-Robertson- Walker (FRW) universes, which provide the stan- 
dard cosmological models, are drastically diluted by the universal expansion, especially during 
inflation. This leads to residual -B-fields that are far too weak to seed and sustain the galactic 
dynamo. For decades, the solution has been usually sought outside Maxwellian electromag- 
netism, standard cosmology, or beyond the linear regime. We can do this in many ways, which 
explains the number and the variety of the proposed scenarios (see [3] for relatively recent work 
and 4J for more references). Nevertheless, one could still produce appreciable magnetic seeds 
of cosmological origin at the linear level and within conventional physics, by appealing to a 
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purely general relativistic (geometrical) amplification mechanism that operates when the FRW 
background universe has negatively curved spatial sections [6]. 



Electromagnetic fields are, so far, the only known energy sources of vector nature. This 
guarantees a purely geometrical coupling between the Maxwell field and the geometry of the 
host spacetime, ensuring a special status for electromagnetism in Einstein's theory. The gravito- 
electromagnetic interaction is monitored by the Ricci identities and adds to the standard inter- 
play between matter and geometry of the field equations. Technically, it ensures the presence 
of curvature-related terms in the electromagnetic wave equations [7]. Physically, this implies a 
very different evolution for the Maxwell field in curved spacetimes. Thus, electromagnetic fields 
in Friedmann models with nonzero spatial curvature evolve differently than in their spatially 
flat counterpart. More specifically, in FRW universes with negative 3-curvature the adiabatic 
magnetic decay slows down and the field is superadiabatically amplified [6]. This is possi- 
ble, despite the conformal invariance of standard electromagnetism and the conformal flatness 
of the FRW spacetimes, because Friedmann models with non-Euclidean spatial geometry are 
only locally conformal to Minkowski space (e.g. see [8]). Thus, the aforementioned amplification 
occurs on relatively large scales, where the curvature effects are prominent and the conformal 
mapping between the Friedmann and the Minkowski spacetimes breaks down. The mechanism 
seems to work in all negatively curved FRW models, irrespective of their matter content and, 
in principle, can lead to residual magnetic fields strong enough to seed the galactic dynamo. 
As yet, this is the only scenario where superadiabatic magnetic amplification is achieved at the 
classical level and without introducing any new physics. 

The purely geometrical mechanism of magnetic amplification outlined above, has been dis- 
cussed and developed within the framework of the Friedmann models in [6] . It applies to more 
general spacetimes, however, since analogous exact results have also been obtained in Bianchi 
class B models that contain electromagnetic fields [9]. Overall, the amplification effect appears 
to be a generic feature of cosmological models with negative spatial curvature. As a result, 
the latter can sustain large-scale magnetic fields much stronger than generally expected. More 
specifically, the final strength of the amplified -B-field can reach magnitudes close to 10~ 15 G 
today, depending on the precise model of inflation [6]. These numerical estimates correspond to 
magnetic modes with sizes close to the curvature scale of an open FRW universe. In particu- 
lar the aforementioned studies followed the evolution of the mode from the moment it crossed 
outside the Hubble radius, during a period of slow-roll inflation, throughout the subsequent 
radiation and dust epochs, all the way to the present. Here, we will turn our attention to the 
early stages of the magnetic evolution prior and during the slow-rolling phase. Our aim will be 
to examine whether the superadiabatically amplified magnetic modes can be causally connected 
at the onset of inflation and the time they crossed the Hubble horizon during the inflationary ex- 
pansion. Causality will ensure the physical existence of these modes, at least in classical (purely 
general relativistic) terms. The time of first horizon crossing, on the other hand, will essentially 
determine the modes' subsequent amplification and eventually their residual strength today. We 
then provide the spectrum of the present amplitude of the superadiabatically amplified magnetic 
modes, both analytically and numerically. According to our results, for reasonable values of the 
cosmological parameters, the final 5-field lies between 10 -20 and 10 -12 Gauss, which means 
that it falls fairly comfortably within the typical galactic-dynamo requirements. 
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2 Magnetic fields on open FRW backgrounds 



Cosmological magnetic fields are thought to decay adiabatically in all FRW backgrounds at 
all times. The reason is believed to be the conformal invariance of classical electromagnetism, 
combined with the conformal flatness of the Friedmann models. Nevertheless, this is not the 
case when the FRW host has hyperbolic spatial geometry. 

2.1 The gravito-electromagnetic interaction 

Electromagnetic fields are, as yet, the only known energy source of vector nature. This means 
that the Maxwell and the gravitational fields have a twofold interaction. The first is the familiar 
interplay between matter and geometry introduced by Einstein's equations. The second is a 
purely geometrical coupling, which holds even outside the realm of general relativity and is 
monitored by the Ricci identities. When applied to the magnetic vector (B a ), the latter read 

2V [a V b] B c = R abcd B d and 2D [a D b] B c = K abcd B d , (1) 

where the former applies to the whole spacetime and the latter holds on the observer's 3- 
dimensional (irrotational) rest-space. Note that V a and D a are the 4-D and 3-D covariant 
derivative operators respectively; R a bcd and Tlabcd are the corresponding Riemann tensors. 

This interaction between gravity and electromagnetism implies that the evolution of large- 
scale cosmological magnetic fields is affected by the geometry of their host spacetime. More 
specifically, the use of the Ricci identities adds curvature-related terms to the wave equation of 
the electromagnetic field [7\. In particular, on an FRW background, the magnetic component of 
a source-free electromagnetic field is monitored by 

B a - B 2 B a = —5HB a -2(H + m 2 + oT 2 K) B a , (2) 

to linear order. Here, H = a/a is the background Hubble parameter (a = a(t) is the associated 
metric expansion scale factor), K = 0, ±1 is the 3-curvature index of the underlying FRW model, 
D 2 = D a D a is the covariant Laplacian operator of the spatial hypersurfaces and overdots denote 
comoving proper-time derivatives. Note the gravito-magnetic term on the right-hand side of 
Eq. ([2]). It results from the magneto-geometrical interactions manifested in the Ricci identities 
(particularly in {lb)) and will play the key role in the analysis that follows. 

We can simplify expression ([2]) by introducing the rescaled magnetic field vector B a = a 2 B a 
and by using conformal, instead of proper, time (i.e. r), with f\ = 1/a). Then, Eq. ([2]), written 
for the n-th magnetic mode, reduces to 

B'( n) + (n 2 + 2K)B {n) =0, (3) 

where n is the eigenvalue of the Laplacian that represents the comoving wavenumber of the 
mode and primes denote conformal-time derivatives [6]. On FRW backgrounds with Euclidean 
spatial geometry (i.e. when K = 0), one can easily show that the rescaled magnetic field remains 
constant in terms of conformal time. This ensures an adiabatic decay for all magnetic modes 
(i.e. Bf n \ oc a -2 , with n > 0) at all times. Following ([3]), for modes with n 2 S> 2, the adiabatic 
decay-rate persists in Friedmann models with nonzero 3-curvature as well. This is to be ex- 
pected, since the aforementioned modes correspond to small scales where the curvature effects 
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are unimportant. On large enough lengths, however, the adiabatic decay law is no longer guaran- 
teed due to the magneto-geometrical term on the right-hand side of ([3]) As mentioned earlier, 
the latter reflects the purely general relativistic interaction between magnetism and spacetime 
curvature monitored by the Ricci identities. 



2.2 Superadiabatic magnetic amplification 

During a period of inflationary expansion the universe is believed to behave like a very poor 
electrical conductor. After inflation and reheating, the conductivity is high and the resulting 
electric currents will freeze any large-scale magnetic field that may be present into the primordial 
plasma. In highly conductive environments the -B-field drops adiabatically at all times (i.e. B a oc 
o -2 ) and irrespective of the background geometry. On scales lying beyond the Hubble radius, 
however, causality ensures the absence of currents. There, the conductivity remains low and the 
ideal magnetohydrodynamic (MHD) approximation does not hold. In such an environment, the 
magnetic field vector is still monitored by Eqs. Q and ([3]). 

According to ([3]), the adiabatic magnetic decay persists on closed FRW backgrounds (where 
K = +1). There, the magneto-curvature term simply modifies the frequency of the field's 
oscillation. On open FRW backgrounds, however, the same term can change the nature of the 
magnetic wave equation. In particular, setting K = — 1 in ([3]), the latter reads 

B'( n) + (n 2 -2)B {n) =0, (4) 

with n 2 > 0. Therefore, on large enough scales (those with —2 < n 2 - 2 < 0), the solution of 
Eq. (j3J) no longer leads to conventional wave solutions, but to 'exponential waves' of the form 

£>( n ) = Ci cosh (rj\/ 2 — n 2 ^ + C2 sinh (j]\/2 — n 2 ^j , (5) 

where now < n 2 < 2. When our open FRW background contains a single perfect fluid of 
barotropic index w = p/p ^= — 1/3, the associated scale factor evolves as a@ oc sinh(/3r/), where 
ft = (1 + 3w)/2. Then, expression ([5]) recasts into the evolution law 



for the actual magnetic field. The above implies superadiabatic amplification, namely a decay 
rate slower than the adiabatic one, for all modes with < n 2 < 2. Also, given that n 2 = 1 
corresponds to the curvature scale of the FRW background (see § 13. II below), we deduce that the 
amplification applies to all supercurvature modes (i.e. those with < n 2 < 1), as well as to the 
largest subcurvature modes (i.e. the ones with 1 < n 2 < 2). More specifically, on the curvature 
scale the magnetic field drops as B oc a , while at the homogeneous limit (i.e. as n 2 — > 0) its 
decay rate slows down further to B oc a/ 2 ^ 2 . 



1 The Minkowski-like evolution of the rescaled B-field on flat Friedmann models can also be seen as the direct 
consequence of the conformal invariance of standard electromagnetism and of the fact that the spatially flat FRW 
spacetime is globally conformal to the Minkowski space. The latter does not apply to Friedmann models with 
nonzero 3-curvature, which are only locally conformal to the Minkowski space [S]. In particular, the conformal 
mapping between the two spacetimes breaks down on large enough scales where the curvature effects become 
prominent. It is on these wavelengths that the magneto-geometrical term in Eq. ((3| takes over and the adiabatic 
magnetic decay law is no longer guaranteed, despite the conformal invariance of standard electromagnetism. 
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So far, we have established that magnetic fields on spatially open FRW backgrounds can 
be superadiabatically amplified by curvature effects within the realm of classical electromag- 
netic theory. We have also identified the scales spanned by the affected magnetic modes and 
found that the amplification effect is essentially independent of the type of matter that fills the 
universe. This means that large-scale -B-fields can be superadiabatically amplified during slow- 
roll inflation, reheating, and throughout the subsequent radiation and dust eras. On physical 
grounds, we would also like to know whether the scales of interest can be causally connected at 
the onset of the inflationary expansion, as well as the time they cross outside the Hubble horizon. 
To large extent, the latter will decide the residual strength of the associated superadiabatically 
amplified -B-fields. It will therefore help to investigate how the hyperbolic spatial geometry of 
the open FRW models affects causality and modifies the standard slow-roll scenario. 



3 Causality and inflation in open FRW models 

When applied to spatially flat Friedmann universes, typical slow-roll inflation leads immediately 
to the familiar de Sitter phase of exponential expansion. This standard picture changes, however, 
in models with non-Euclidean spatial geometry, especially during their early stages. In what 
follows, we will discuss certain aspects of FRW cosmologies with negative 3-curvature. 



3.1 Causal horizons 

Consider a Friedmann universe with hyperbolic spatial hypersurfaces, that contains a single 
barotropic medium of energy density p and isotropic pressure p. The evolution of this model is 
determined by the set of equations 

p = -3flp(l + w) , H = -H 2 - - k P (1 + 3w) (7) 



and 



r2 1 1 



H z = - K p + -, (8) 

o (X 

where w = p/p is the barotropic index and k = 8nG is the rescaled gravitational constant (with 
c = 1). Recalling that Xh = 1/H is the Hubble horizon, Q = np/3H 2 is the density parameter 
and defining \k = a as the curvature scale of the universe, we see that Eq. ([8]) is recast as 

with < O < 1 for all open FRW models (e.g. see [IQ]). The latter means that the curvature 
radius of an open FRW cosmology lies always outside the Hubble radius, with the two scales 
getting close in the low density limit (i.e. \k — > A^- as Q — > + ). 

Let us now go back to the magnetic case discussed in § 12.21 and consider a mode with physical 
wavelength A ra = a/n. We saw that large-scale -B-fields are superadiabatically amplified during 
inflation, reheating, and throughout the subsequent evolution of the universe. Given the large 
wavelengths of these modes, it is of interest to establish whether or not they can be causally 
connected at the onset of inflation. Recalling that \k = a, and using @, we arrive at 

A - = Af (10> 
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which relates the scale of the magnetic perturbation to that of the Hubble horizon. The above 
ensures that the scale of the mode in question will lie inside the Hubble radius, and therefore it 
will be in causal contact, as long as 

n<i-\. (ii) 

This condition is readily satisfied by all subcurvature lengths, namely by those with n 2 > 
1. In other words, scales smaller than the curvature radius of a open FRW spacetime can 
reside inside the Hubble horizon, provided the universe is sufficiently open. These include the 
largest subcurvature magnetic modes (i.e. those with 2 > n 2 > 1), which are superadiabatically 
amplified (see § I2.2[) . For example, subcurvature magnetic modes with n = 1.01 and n = 1.1 
will be within the Hubble radius if SI < 0.02 and Q < 0.17, respectively, at the beginning of the 
inflationary phase. Supercurvature modes can never satisfy condition (jlip . This reflects the fact 
that the corresponding lengths are always outside the Hubble radius. Nevertheless, in principle, 
even supercurvature lengths can be causally connected in open FRW models. 

When dealing with spatially flat FRW cosmologies containing conventional matter, namely 
with w > —1/3, the Hubble length essentially coincides with the particle horizon of the comoving 
observers and therefore it effectively defines their observable universe. In those models, super- 
Hubble (and consequently super-curvature) scales are always causally disconnected. This is not 
the case, however, in open Friedmann models. There, the Hubble and the particle horizons do 
not coincide, but the latter is generally larger than the former. For instance, at any given time, 
the particle horizons of a radiation and a dust-dominated open FRW universe are 



A f= ' h l±P and A„ = -i=ln (1 + ^3 , (12) 

respectively (see also [II]). Then, a simple Taylor expansion shows that both particle horizons 
diverge as the energy density of the universe becomes progressively lower (i.e. Xp — > +00 as 
£1 — > + in agreement with the Milne limit). Close to the Euclidean threshold, on the other 
hand, one recovers the familiar Xp = 1/H and Xp = 2/H expressions of the radiation and the 
dust eras respectively. The conclusion is that super-Hubble, as well as super-curvature, scales 
can be causally connected in FRW universes with sufficiently low densityH It is therefore of 
physical interest to explore their phenomenological consequences. 

3.2 The inflationary dynamics 

Introducing conformal time and assuming a period of de Sitter- type inflation (when p = —p) 
the system of ([7]) and ([8]) has the parametric solution 



a = ao 



e V(l _ e 2r?o) 

e vo(i - e 2r i) 



(13) 



2 Our approach is purely general relativistic and causality is the criterion that decides whether physical processes 
can operate, or not, within a given region. It has been argued that, quantum mechanically, the superadiabatically 
amplified modes may not be normalisable 12 . The argument was based on a study adopting open inflation, instead 
of the standard slow-roll scenario used here. The authors found that the modes were normalisable on small scales, 
where the 3-curvature effects are negligible, but encountered the usual infinities when they reached lengths where 
the curvature is strong. Leaving aside the question of open inflation, the aforementioned normalisability issues 
are not surprising, given the absence of a theory unifying general relativity and quantum mechanics. After all, 
these are exactly the problems that any future theory of quantum gravity is expected to solve. 
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and 



u 2e»w 



(l + e^l-e^ 



(14) 



.(1 - e r ')(l + evo) 

where r\ < and the zero suffix denotes a given time during inflation. Thus, at the beginning of 
the inflationary phase, when r\ <C and e'<l, the scale factor evolves as a oc t. Towards the 
end of inflation, on the other hand, we have r/ — > CP and — ¥ 1", which ensure that a oc e*. Put 
another way, the slow-rolling regime in an open FRW model starts with coasting expansion and 
the exponential phase only occurs at the end, when the effects of curvature have faded away. If 
we are interested in modes close to the curvature scale, it is important to quantify the effect of 
the coasting phase. 

Further insight can be obtained by evaluating the principal kinematic and dynamic variables. 
In particular, using primes to denote conformal time derivatives, the Hubble and the deceleration 
parameters are H = a' /a 2 and q = — (l + H'/aH 2 ) respectively. Then, on using expressions (fl~3|) 
and (I14p . we obtain 

H= a eV(l-e 2 vo) > &nd 9 = "(TT^F' (15) 
The former combines with Eqs. (jHJ) and (I13p to give 

1 4e 2r)0 

3 KP= ag(l-e*»)*' (16) 

which guarantees that p = po = constant, as expected. Finally, substituting this result in the 
right-hand side of (|15b) and keeping in mind that Q = Kp/3H 2 , we arrive at 

2 V 3 y v 1 

Accordingly, at the onset of inflation (i.e. when r) —> — oo) we have H — > +oo, q — > _ and 
O — > + . On the other hand, towards the end of the inflationary expansion (i.e. as rj —> 0~), 
we find that H — > y / Kpo/3, O — >• 1 _ and g — )• — 1 + . As before, we see that the model enters its 
accelerating, de Sitter-type, phase only asymptotically at the end of the inflationary regime. 



3.3 Hubble horizon crossing 

During the coasting phase of the expansion, we have H oc t , which means that the Hubble 
horizon scales like Xh oc a. Given that all wavelengths also scale in the same way, modes that 
were originally inside the Hubble radius will remain so and cross outside only towards the end 
of inflation, when the curvature effects have faded away and the expansion starts to accelerate. 
As will see next, the time of horizon crossing is decided solely by the comoving wavelength of 
the mode in question. 

Let us consider a mode with physical scale A n = a/n, where n is the associated comoving 
eigenvalue. On using expressions (fl~3j) and (fl~5h). we obtain 

\ L = aH = 1 + e 2 ^ 
Xh n n(l — e 2r ?) ' 
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with rj < 0. Thus, a mode will lie inside the Hubble radius (i.e. X„/Xh < 1) as long as 

1 + e 2r > 

n> T^' (19) 

Given that r/ < 0, the latter guarantees that n > 1 always. In particular, the right-hand side of 
the above approaches unity at the beginning of the inflationary phase (when r/ — > — oo) and starts 
to diverge as we approach the end of inflation (i.e. as r] — > 0~). In other words, supercurvature 
scales (those with < n 2 < 1) are always outside the Hubble horizon. 

Taking the conformal time derivative of (|18p . we find that (A n /A//)' = 4e 2?? /n(l — e 2?? ) 2 > 
at all times, which ensures that during the inflationary regime any given wavelength grows 
faster than the Hubble horizon. This in turn implies that, given enough time, essentially all 
sub-Hubble scales will eventually cross outside. Following (|19p . at the time of horizon-crossing, 
namely when A n /A# = 1 and rj = tjhc, we find that 

1 , fn- 1\ 

tjhc = - In — — < , (20) 



2 \n+l, 

since n > 1 (see Eq. (|19p above). Consequently, all finite wavelengths can exit the Hubble scale 
before the end of inflation. However, for modes that are initially well inside the Hubble radius 
(i.e. with n> 1), horizon crossing occurs only at the end of inflation (i.e. t\hc ~ > 0~ ). Further 
information can be obtained by recasting the Friedmann equation into the form 

Xh n\/l — Q, ' ^ 

with < Q < 1 in our case. Therefore, at the time when a given wavelength passes through the 
Hubble threshold, the density parameter of the universe is 

Sl HC = 1 - \ ■ (22) 

where n 2 > 1@ A gain, we see that most subhorizon modes do not cross the Hubble scale before 
the value of the density parameter has approached substantially close to unity. Wavelengths 
spanning scales just inside the curvature radius, with n = 1.01 and n = 1.1 for example, cross 
outside the Hubble horizon at Qhc = 0.02 and Qhc = 0.17 respectively. For smaller lengths, 
say with n = 2 and n = 10, horizon crossing occurs later, when £Ihc = 0.75 and flsc = 0.99 
respectively. Finally, it is worth pointing out that a given mode re-enters the Hubble horizon 
(during the subsequent radiation or matter dominated eras) when Q is again equal to £Ihc- For 
the current Hubble scale, whatever n it possesses, this means that Qhc = a result that will 
prove useful later. 

3.4 The number of e-folds 

In a FRW universe, the Hubble scale is related to the density parameter by means of the 
expression Xh = y/3£l/K,p. Therefore, during inflation, a perturbation that crosses the Hubble 
horizon has wavelength 



if 

M 2 



X n = X H ^ — , (23) 



3 According to expression (|22f) . modes with < n 2 < 1 have Q,hc < 0. This result simply reflects the fact 
that the corresponding scales lie outside the Hubble length at all times. Recall that n 2 = 1 corresponds to the 
curvature radius, which in open FRW models is always larger than the Hubble horizon (see § [371] earlier) . 
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with n 2 > 1. Note that k ~ Mp? and p ~ M 4 = constant, where Mpi represents the Planck mass 
and M the energy scale of inflation in natural units. After horizon crossing, the aforementioned 
wavelength grows by a factor of e , where N is the number of e- folds between horizon-crossing 
and the end of the inflationary regime, so then we have 

(\ \ „ „N M Pl n HC (0A \ 

Throughout the reheating era, the energy density of the matter decreases as p oc a~ 3 . Conse- 
quently, by the end of that period, the scale in question has grown further by 



a RH _ fpiNF\ l/Z _ ( M \ 4 / 3 
\ PRH ) \TrhJ 



(25) 



aiNF 

since prh — T RH (in natural units). Putting (|23j) and ff2fi§ together, we deduce that at the 
beginning of the radiation era, the original wavelength is 

(Xn)nH^e N ^%, (26) 

with the quantities in the right-hand side given in natural units. 

Expression (|26p provides the scale spanned at the end of reheating by a wavelength that 
crossed the horizon N e-folds before inflation was over. We may also estimate the aforementioned 
number of e-folds using thermodynamic arguments. More specifically, following [13j . we assume 
that the expansion of the universe proceeds adiabatically (with the exception of the reheating 
era). Then, the entropy contained within the current Hubble scale, which is estimated close to 
10 88 , has remained unchanged since reheating. At that time, the entropy inside a region that 
crossed the horizon N e-folds before inflation ended was 

RH - x RhTrh - e M 2 Trh . ( 27 ) 



where the last equality derives from Eq. (|26]h Putting Srh equal to 10 in the above and 
recalling that M Pi ~ 10 19 GeV, we find that a scale of the size of the current observable 
universe, crossed the Hubble horizon 

AT* ~ 24 + | InM + i lnT^-i mfi* H C, (28) 

e-folds before the end of inflation, where the * denotes quantities associated with the current 
Hubble scale of the universe. 

In general, a wavelength of size X n today crossed the Hubble horizon N e-folds before the 
end of inflation when rj = r?#c, with the latter given in Eq. (|20|) . Obviously, wavelengths smaller 
than the present Hubble scale crossed later and larger scales crossed earlier. Then, using (fT3j) , 
we have 

n* 1 — P 2 VHC 

e N-N* = ^HC = e v' BC -VHO 1 6 . (29) 

a,HC 1 — e v hc 

Finally, substituting (1281 into the above, while employing ([20]) and (f22j) . yields 



jV " 53 + i'"(i^) + 5 ln (i^)-5 ln < 1 - !! » ) -5 ln (" 2 - 1 )- ( 3 »> 
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given that = fin (see Eq. (|22p in § 13. 3p . In the above, which generalises Eq. (8.45) of [13] 
for FRW models with Q < 1, the quantities M and Trh are measured in GeV. Also,n is related 
to the physical wavelength at present by n = (1 — ^o) _1//2 (AffMn)o ( see a ls° (|^T]) )n 

Looking at expression (|30f) . one can see that the number of e- folds does not scale as N oc In X n , 
in contrast to the flat-FRW case (compare to Eq. (8.45) in [13]). The difference reflects the fact 
that, in a open FRW models, the Hubble radius is not constant during inflation (see § 13.21 
earlier). Also, following (|30p . the e- folding number diverges as the wavelength approaches the 
curvature scale (i.e. for n 2 — > 1). However, the singularity is apparent and can be explained 
by recalling that the supercurvature scales, namely those with < n 2 < 1, are always outside 
the Hubble length. Hence, for the corresponding modes, there is no horizon exit and Eq. (|30p 
does not apply. In their case, the key parameter is the total number of e-folds. Estimating the 
amplification of supercurvature S-modes is not essential for our purposes and goes beyond the 
scope of the present article. Therefore, from now on, we will focus on magnetic fields spanning 
the largest subcurvature scales (i.e. modes with 1 < n 2 < 2). 



4 The residual magnetic field 

The number of e-folds between the time a certain magnetic mode crossed outside the Hubble 
radius and the end of inflation, is crucial for the current strength of the residual i?-field. There 
are additional factors, however, that can also play a key role. 



4.1 The overall amplification 

In line with the Gibbons-Hawking temperature, the energy density (pb) stored inside a magnetic 
mode at horizon crossing, is determined by the dimensionless ratio 



Pb\ (MY 
\M Pl J 



P J HC 



(31) 



P / INF 



where pb ~ B 2 is the energy density of the S-field. Magnetic fields coherent on the largest 
subcurvature scales (and beyond) are superadiabatically amplified as B oc a~^ n \ with (3(n) = 
2 — \/2 — n 2 . Then, given the constant background density throughout inflation, we have 

■,B\ ~ e -«»>(^)\ (32) 

at the end of the slow-rolling regime. 

During reheating, the effective equation of state of the matter is that of pressure- free dust. 
This means that pre = Pinf{o-inf / clrhY by the end of that period. At the same time, the large 
scale magnetic fields are still superadiabatically amplified. Therefore, near the curvature scale, 
where B oc a~P( n \ we have (pb)rh = {pb)inf{o-inf I '(Irh) 2 ^^ '■ As a result, since pinf — M A 
and prh — Tr H , the relative strength of the associated magnetic field at the end of reheating is 




RH 



\Trh 



) [wJ ■ (33> 



4 Substituting the relation n = (1 — fio) 1 ^ 2 (Af//A n )o to the right-hand side of Eq. (|30[) and then setting Qq = 1, 
one can easily recover expression (8.45) in [13] . 
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Recall that prh — P-y at that time, with p 1 representing the energy density of the radiative 
component. The magnetic amplification continues throughout the radiation epoch and the 
subsequent dust era. Therefore, for 5-fields spanning lengths close to the present curvature 
scale of the universe, we have (pb/p^o = {pb/ ' P~i)rh(Trh /Tn) 2 [ 2_/3 ( n ^ today. Also, to obtain 
the last result, we have taken into account that p 7 oc a -4 and T oc 1/a. Combining these with 
(1301) and setting T ~ 10" 13 GeV, gives 




-LQ-117+102V2^n 5 ' 



Jqu) [(l-^o)(n 2 -l)] , (34) 



since /3(n) = 2 — y/2 — n 2 and Vt* HC = f2 (see Eq. (f22j) in § I3.3[) . Finally, recalling that 
(py)o ~ 10~ 51 GeV 4 , we arrive at an expression for the present-day magnetic field 



B Q ~ W -G5+5lV2^ 



M 



2VT- 



„ 1 (2-v / 2-n 2 )/2 

l-n )(n 2 -l) V , (35) 



TO 14 , 

for the largest subcurvature modes (i.e. those with 1 < n 2 < 2). This formula provides, in Gauss 
units, the present magnitude of a superadiabatically magnetic mode that spans scales close to 
the curvature radius of an open FRW universe. 



4.2 The final magnetic spectrum 

According to expression ([35]) . the final strength of the superadiabatically magnetic field depends 
on the energy scale of the adopted inflationary scenario (M), on the current density of the 
universe (On) and on the mode's scale (n). The WMAP-normalised data constrain the current 
density parameter of the universe, so that approximately 1— Q,q < 10~ 2 . Also, typical inflationary 
models have energy scales between 10 14 GeV and 10 17 GeV. Thus, setting 1 - fi ~ 10" 2 and 
M ~ 10 14 GeV, the current magnitude of a (subcurvature) magnetic mode with n = 1.01 is 

B ~ 10~ 16 Gauss. (36) 

Stronger i?-fields can be obtained by increasing the scale of inflation. For instance, keeping 
the wavelength of the mode and the current density parameter unchanged and assuming that 
M ~ 10 16 GeV, we arrive at -Bo ~ 10~ 12 G. Magnetic seeds of such strengths are still too 
weak to affect outcome of primordial nucleosynthesis [13], or leave an observable signature in 
the CMB [15J. They can readily sustain the galactic dynamo, however, which typically requires 
an initial B-seed with magnitude in the range between 10 -20 G and 10 -12 G [2]. In addition, 
intriguingly, magnetic fields close to 10~ 16 G are very close in strength to those recently reported 
in empty intergalactic space [3]H 

The numerical estimates quoted above will drop if the current density parameter gets closer 
to unity. However, the dependence of the residual magnetic strength on the present density 

5 When all the magnetic amplification takes place during inflation and the B-field decays adiabatically ever 
since, strengths around 10 -16 G today can lead to the so called backreaction problems [16]. These occur because 
the required inflationary amplification is so strong that the magnetic energy density catches up with that of the 
driving inflaton. In our scenario there is no backreaction issue. The energy density of the B-field is much lower 
than that of the dominant matter component at all times, namely during inflation, reheating and later during 
the radiation and the dust epochs. Residual strengths close to 10~ 16 G are achieved because the superadiabatic 
amplification of the field persists throughout its evolution and it is not confined to the inflationary era only. 
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Figure 1: The current spectrum of the superadiabatically amplified magnetic field, as a function 
of the mode wavelength, for different energy scales of inflation, and values of the current density 
parameter (on a logarithmic scale). Small wavelengths experience weak (or no) amplification, 
while the maximum enhancement occurs very close to the curvature scale. The cut-off seen 
at the curvature radius of the universe simply reflects the fact that supercurvature modes are 
not included in our analysis. Note that, as Qq approaches unity, the curvature scale is pushed 
further away from the Hubble horizon, which is assumed to be at approximately 10 3 Mpc. The 
dependence of the numerical results on Qq, however, is considerably weaker than that on M. 

of the universe is relatively weak, and certainly weaker than that on the energy scale of the 
adopted inflationary model. This means that we can obtain astrophysically relevant 5-seeds 
even in very marginally open universes. For example, setting n = 1.01, M ~ 10 14 GeV and 
1 - O ~ 10~ 10 , we find B ~ 10" 20 G. For M ~ 10 16 GeV, the same strength can be obtained 
even when 1 — ~ 10~ 18 . These results correspond to i?-modes coherent on lengths just 
inside the curvature radius of an open FRW universe, which experience strong superadiabatic 
amplification. In particular, a mode with n = 1.01 decays as B oc a -101 (see solution ([6]) in 
§ I2.2p . As we move down to smaller scales, however, the magnetic decay-rate increases and 
at the n = 2 threshold the adiabatic decay law is re-established (see Fig. [I]). We also remind 
the reader that supercurvature lengths, although causally connected, are not included in our 
analysis (see § 13. 4p . 

It should be emphasised that the above quoted magnetic strengths do not include the am- 
plification of the field during the nonlinear phase of galaxy formation. In the case of spherically 
symmetric protogalactic collapse, for example, the magnitude of the magnetic seed typically 
increases by two or three orders of magnitude. In the more realistic case of anisotropic collapse, 
shearing effects could add one or two extra orders of magnitude to the .B-seed [17]. Then, the 
galactic dynamo is expected to take over [2]. Also, our numerical results refer to S-fields span- 
ning scales comparable to the curvature radius of the universe. Once galaxy formation begins, 
however, these fields should start breaking up and reconnecting on scales comparable to that of 
the collapsing protogalaxy. Note that, on sufficiently small scales, damping effects will process 
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the magnetic spectrum and determine its small-scale structure and polarization |18j . 

5 Discussion 

The problem of cosmic magnetism - its origin, evolution and implications - persists. The con- 
tinual detection of magnetic fields in the universe, makes the search for its solution increasingly 
pressing. For cosmologists, the biggest challenge is obtaining sustainable primordial -B-fields 
that could successfully seed galactic dynamos. There have been many attempts to solve this 
problem. With very few exceptions, however, the solution has been sought outside standard 
electromagnetism, standard cosmology, or general relativity. Nevertheless, the resources of Ein- 
stein's theory may not have been exhausted yet and it could still provide the answer. The 
reason is the geometrical interpretation of general relativity, which gives a special status to the 
electromagnetic field, due to its vector nature. In practice, this means that, in addition to the 
Einstein equations, the Maxwell field interacts with the geometry of the host spacetime through 
the Ricci identities as well. Although this particular coupling between gravity and electromag- 
netism has been known for some time, it remains under-investigated. Here, we have explored the 
implications of the aforementioned interaction for the evolution and the survival of large-scale 
cosmological magnetic fields in greater detail. 

We have studied how the generation of seed magnetic fields during standard slow-roll inflation 
in a marginally open universe can lead to significant superamplification of those seed fields to 
levels that can provide the initial fields needed for dynamos to produce easily the -B-fields with 
magnitudes now observed in galaxies, clusters and intergalactic space. This process of magnetic 
superamplification relies crucially on the existence of negative spatial curvature so that the 
primordial magnetic fields do not decay away rapidly at the adiabatic, B oc a~ 2 , rate. Our 
mechanism appeals to the aforementioned magneto-curvature coupling, which in open universes 
can significantly slow the decay of the magnetic field, by breaking the global conformal invariance 
that has often wrongly been assumed to dominate its evolution. The resulting superadiabatic 
magnetic amplification applies to .B-fields coherent near and beyond the curvature radius of the 
universe. Here, we have looked at the largest subcurvature scales, calculated the fields expected 
for likely parameter choices and displayed the predicted magnetic spectrum (see Fig. [I]) . In all 
cases, the amplification peaks just inside the curvature scale, which implies a typical correlation 
length of the order of 10 4 Mpc for the maximally amplified -B-field. The latter, however, should 
break up and reconnect on much smaller (cluster or galactic-size) scales once the nonlinear phase 
of galaxy formation begins. The residual magnetic strength depends on the energy scale of the 
adopted inflationary model and on the current curvature radius of the universe. More specifically, 
setting the energy scale of inflation at M ~ 10 14 GeV 4 and assuming that 1 — Qo ~ 10 -2 , the 
maximum magnetic strength is approximately 10~ 16 G. 

These numerical results can be tuned further when the data from the Planck mission narrow 
down the range of Q,q and that of the inflation energy scale. Until then, we should underline that 
in our scenario the earlier inflation starts, the stronger the magnetic amplification. On the other 
hand, the closer the density parameter gets to unity, the weaker the final .B-field. Nevertheless, 
the ^-dependence is relatively weak, which means that even in (very) marginally open universes 
the residual magnetic seeds are comfortably within the standard dynamo requirements. Finally, 
before closing, we should also point out that -B-fields around 10 -16 G are intriguingly close to 
those recently reported in empty intergalactic space. Overall, our analysis shows that spacetime 
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curvature can act as an effective dynamo mechanism of purely geometrical nature. 
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